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Abstract: In general the gravitino mass and/or the soft supersymmetry breaking 
masses in the observable sector can be much larger than the TeV scale. Depending 
on the relation between the masses, new important channels for gravitino production 
in the early Universe can arise. Gravitinos with a mass above 50 TeV decay before 
big bang nucleosynthesis, which leads to relaxation of the well known bound on the 
reheating temperature Tr < 10^° GeV. However, if the heavy gravitinos are produced 
abundantly in the early Universe, their decay can alter the abundance of the lightest 
supersymmetric particle. Moreover, they may dominate the energy density of the 
Universe. Their decay will in this case increase entropy and dilute already created 
baryon asymmetry and dark matter. Such considerations put new constraints on 
gravitino and sfermion masses, and the reheating temperature. In this paper we 
examine various cosmological consequences of supermassive gravitinos. We discuss 
advnatges and disadvantages of a large reheating temperature in connection with 
thermal leptogenesis, and find that large parts of the parameter space are opened 
up for the hghtest right-handed (s)neutrino mass. We also discuss the viabihty of 
Affleck- Dine baryogenesis under the constraints from gravitino decay, and gravitino 
production from the decay of Q-balls. 
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1. Introduction 

Primordial inflation is the most convincing paradigm for the early Universe |l| . The 
vacuum fluctuations created during inflation also explain the observed temperature 
anisotropy of the cosmic microwave background (CMB) radiation p|. However in- 
flation leaves the Universe cold and void of any thermal entropy. Entropy is be- 
lieved to be created from the decay of the coherent oscillations of the inflaton which 
can happen perturbatively and/or non-perturbatively into bosons ^, ^ and 
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fermions [0]. It is necessary that the standard model (SM) degrees of freedom are 
produced at this reheating stage, particularly baryons which are required for the 
synthesis of light elements during big bang nucleosynthesis (BBN) at a temperature 
O(MeV) 



On the other hand, we do not know the full particle content of the Universe be- 
yond the electroweak scale, therefore we do not know what degrees of freedom were 
excited right after inflation. In this respect, supersymmetry (SUSY) acts as a build- 
ing block which can explain the hierarchy between the Planck and electroweak scales. 



if it is softly broken in the observable sector at the TeV scale, see [0. Besides its 
phenomenological implications, this also has important cosmological consequences. 
The scalar potential of the minimal supersymmetric standard model (MSSM) has 
nearly 300 flat directions [0. These fiat directions can address cosmological issues 



from reheating to density perturbations [T^. If i?-parity is conserved, the lightest su- 
persymmetric particle (LSP) will be stable and can act as a cold dark matter (CDM) 
candidate. A neutralino LSP with a mass ~ 100 GeV can match the current 



observational limit when produced thermally |14 



SUSY breaking at the TeV scale in the observable sector can be achieved via 
gravity |TT|, gauge |TB| and anomaly |]TB| mediation, leading to different patterns of 
supersymmetric particle masses. However, there is a priori no fundamental reason 
why this scale should be favored by nature ^. 

Inspired by the string landscape |]T7|, |TB|, there has recently been an interest- 
ing proposal for SUSY breaking well above the electroweak (but below the Planck) 
scale pO| , In this new scheme, coined split SUSY, the masses of sfermions can 
be arbitrarily larger than those of fermions. Although such a scheme does not at- 
tempt to address the hierarchy problem, it removes fear from flavor changing and 
CP— violating effects induced by the light scalars at one-loop level. Successful gauge 
coupling unification requires that the gauginos be kept lighter than 100 TeV, while 
spontaneous breaking of the electroweak symmetry requires that the mass of the 
hghtest Higgs be around C(IOO) GeV. 

A priori there is no fundamental theory which fixes the scale of SUSY breaking, 
but cosmological considerations can constrain it. For example, the theory permits a 



^String theory which is believed to be the most fundamental theory does not provide us with 
a concrete answer. Rather it provides us with a landscape with multiple vacua where the 
SUSY breaking scale remains undetermined p8[ |, this transcends into an uncertainty into the scale 
of inflation and the required number of minimal c- foldings . 
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very long-lived gluino. The annihilation of gluinos alone may not efficiently reduce 
their abundance below the experimental limits on the anomalous isotopes of ordinary 
matter |2y]. The decay of gluinos within the lifetime of the Universe solves this 
potential problem, and requires sfermions masses to be less than 10^^ GeV pT| ]. 
Models which give rise to such a split pattern of SUSY breaking masses are typically 
more complicated than the conventional ones for TeV scale SUSY breaking ^ 



22[| . Of course, one might also give upon supersymmetric gauge coupling unification 
and allow all supersymmetric particles to be much heavier than TeV. In this case 
SUSY will be irrelevant for physics at the electroweak scale. On the other hand, 
even in the context of MSSM, some of the sfermions can have a mass ^ 1 TeV, as 



happens in the so-called inverted hierarchy models Therefore, under general 

circumstances, it is possible that at least some of the sfermions are much heavier 
than TeV. 

Local SUSY naturally embeds gravity, hence supergravity, and implies the ex- 
istence of a new particle known as the gravitino which is the superpartner of the 
graviton which is fairly long-lived. Massive gravitinos consist of helicity ±1/2 (longi- 
tudinal) and helicity ±3/2 (transverse) components. In the early Universe gravitinos 
can be produced from thermal scatterings of gauge and gaugino quanta |2^, ^], 



and from the decay of sfermions [25]. Gravitinos are also produced non-thermally 



from the direct decay of the inflaton |26, 27, Eq], and from the vacuum fluctuations 



during the coherent oscillations of the inflaton field |oO|, |oT], In the minimal 
supergravity models, the gravitino mass ms/2 is the same as the soft breaking mass 



of scalars [11]. However, gravitinos can be much heavier once one goes beyond the 



minimal supergravity, for example in no-scale models ]|3^. It is therefore possible 



that m3/2 ^ 1 TeV, even if SUSY is broken at TeV scale in the observable sector. 

Gravitinos with a mass below 50 TeV decay during and/or after BBN. Depending 
on the nature of the decay, there exist tight bounds on the reheating temperature, 
see [|3^, Obviously these bounds do not apply to supermassive gravitinos, i.e. 
when m3/2 > 50 TeV. In such a case the reheating temperature could potentially 
be as large as the inflaton mass, leading to many interesting consequences. As an 



example, it opens up new regions of the parameter space for thermal leptogenesis [pq] , 
a scenario which is sensitive to the reheating temperature as it requires the excitation 
of the lightest right-handed (RH) neutrinos and their supersymmetric partners from 
the thermal bath. 
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However other considerations can constrain the abundance of supermassive grav- 
itinos. Every gravitino produces one LSP upon its decay. If the gravitino decays after 
the thermal freeze-out of LSPs, then it can alter the LSP abundance. In addition, 
gravitinos can dominate the energy density of the Universe if they are produced 
abundantly. Entropy release from gravitino decay in this case dilutes any generated 
baryon asymmetry. 

However, this may turn into a virtue, for example, in the case of Affleck-Dine 
(AD) baryogenesis ^7^. Depending on the parameter space, the AD mechanism 
(which utilizes supersymmetric flat directions) can generate order one baryon asym- 
metry. This would then be diluted to the observed value if the gravitino decay gener- 
ates enough entropy. Often the flat directions also fragment to form non-topological 
solitons such as supersymmetric Q-balls. The Q-balls decay slowly through their 



surface [^, and can themselves be a major source of late gravitino production. 

In this article we consider various cosmological consequences of models with su- 
perheavy gravitinos and/or sfermions, without delving into model-building issues. 
Such particles will be inaccessible at future colliders, and hence cosmology will be 
essentially the only window to probe and/or constrain these models. The rest of 
this paper is organized as follows. In Section II we briefly review various sources of 
gravitino production and their individual contributions. We then discuss gravitino 
decay and constraints from the LSP dark matter on the abundance of supermassive 
gravitinos in Section HI. We briefly review thermal leptogenesis and the effects of 
the gravitino in Section IV, and identify regions of the parameter space which al- 
low successful leptogenesis. Section V is devoted to supersymmetric flat directions 
baryogenesis, including thermal effects and viability of the AD mechanism in the 
presence of supermassive gravitinos. We also show that long-lived Q-balls can be a 
source for copious production of gravitinos. We summarize our results and conclude 
the paper in the flnal Section VI. 

2. Gravitino Production 

The most important interaction terms of the gravitino fleld ipf^ come from its coupling 
to the supercurrent. In the four- component notation, and in flat space-time, these 
terms are written as (see for instance PU[ ) 
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Here X and x denote the scalar and fermionic components of a general chiral super- 
field, respectively, while Pju) and A*-"^ denote the gauge and gaugino field components 
of a given vector superfield respectively. 

In the limit of unbroken SUSY gravitinos are massless and the physical degrees 
of freedom consist of the helicity ±3/2 (transverse) components. After spontaneous 
SUSY breaking gravitino eats the Goldstino and obtains a mass m3/2 through the 
super Higgs mechanism, and helicity ±1/2 (longitudinal) states appear as physical 
degrees of freedom. When the value of is much smaller than the momentum 
of the gravitino, the wave-function of the helicity ±1/2 components of the gravitino 
can be written as 

^ 1 



2 



3m3/2 



:2.2) 



with if) being the Goldstino. The helicity ±1/2 states of the gravitino will in this 
case essentially interact like the Goldstino and the relevant couplings are given by 
an effective Lagrangian 



eff 



1 ±75 



± h.c. 



(2.3) 

Here mx and denote the mass of X and x fields, respectively, while m\ is the 
mass of gaugino field A*^*^-*. The interactions of helicity ±1/2 and helicity ±3/2 states 
of the gravitino have essentially the same strength when \mx — vay\ and m\ are 
smaller than m3/2. In the opposite limit, the rate for interactions of helicity ±1/2 
states with X and x, respectively gauge and gaugino fields, will be enhanced by a 
factor of [vn^x ~ "^xf I ^'^'^zi2 (-^ being the typical energy involved in the relevant 
process), respectively m^^jm}^!^^ compared to that of helicity ±3/2 states. 

Gravitinos are produced through various processes in the early Universe, where 
the relevant couplings are given in (|2.1|) and (|2? 



Thermal scatterings: 

Scatterings of gauge and gaugino quanta in the primordial thermal bath is 
an important source of gravitino production, leading to (up to logarithmic 
corrections) [p^, [25[]: 



Helicity ± - 



Helicity ± 



^3/2 \ 
S ) i 



1± 



Tt 



R 



12"^3/2, 



IQio GeV 



228.75 



3/2 



10" 



-12 



n 



3/2 



R 



IQio GeV 



228.75 
^*(^r). 



3/2 



10 



-12 



(2.4) 
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where Tr denotes the reheating temperature of the Universe, Mg is the gluino 
mass and 5'*(Tr) is the number of relativistic degrees of freedom in the thermal 
bath at temperature Tr. Note that for Mg < 1713/2 both states have essen- 
tially the same abundance, while for Mg ^ ^3/2 production of helicity ±1/2 
states is enhanced due to their Goldstino nature. The linear dependence of 
the gravitino abundance on Tr can be understood qualitatively. Due to the 
Mp-suppressed couplings of the gravitino, see and ( p.3| ), the cross-section 
for gravitino production is oc Mp^. The production rate at temperature T 
and the abundance of gravitinos produced within one Hubble time will then 
be oc and oc T respectively. We remind that the Hubble expansion rate at 



temperature T is given by if = ^(c/,7r790)TVMp, where is the number 
of relativistic degrees of freedom in the thermal bath at temperature T. This 
implies that gravitino production from scatterings is most efficient at the high- 
est temperature of the radiation-dominated phase of the Universe, i.e. when 



Decaying sfermions: 

If m3/2 < fh, the decay channel sfermion — > fermion+ gravitino is kinematically 
open. When m > few x 7713/2, the decay rate has a simple form 



Helicity ± - : 



Helicity ± - : 



sferm^ferm+ip 



sfervn^ferm+ip 



487r ml^^Ml ' 



487rM2 



:2.5) 



Sfermions will reach thermal equilibrium abundances, provided that Tr > m. 
They promptly decay through their gauge interactions when the temperature 
drops below rh. Gravitinos are produced from sfermion decays for the whole 
duration sfermions exist in the thermal bath t ~ M-p/rh?. The abundance of 



gravitinos thus produced will then be ||25 
1 



Helicity ± 
Helicity ± 



^3/2 \ 
s ) 

^3/2 \ 
s ) 



dec 



dec 



m 

"^3/2 

m 
1 TeV 



m 



1 TeV 
228.75 



228.75 



3/2 



N 
46 



1.2 X 10 



-19 



-1 3/2 



N 
46 



1.2 X 10"^^ 



(2.6) 



where g*{m) is the number of relativistic degrees of freedom at T = m, and 
is the number of all sfermions such that 7713/2 < fh < Tr. This result is 



- 6 - 



independent of Tr, so long as Tr > m. Note that for m > m^i2 gravitinos of 
both hehcities will be produced with approximately the same abundance. 

If m ^ 'm^l2-, helicity ±1/2 states interact with sfermions and fermions very effi- 
ciently and can actually reach thermal equilibrium, thus leading to (?^3/2/s) — 

10~^. This will happen when m > ( iC^mg^j^P ) 5 example, m > 10^ GeV 
if m3/2 ~ 50 TeV. 



Inflaton decay: 

Reheating of the Universe also leads to gravitino production p9| , |30| , |26| , |27[| (for 
related studies, see [B3, HO, ET|). Here we consider the case where inflaton decays 



perturbatively and a radiation-dominated Universe is established immediately 
after the completion of its decay.^ This in general provides a valid description 
of the last stage of inflaton decay, regardless of how fast and explosive the first 
stage of reheating might be due to various non-perturbative effects M, H. 



Let us denote the SUSY-conserving mass of the infiaton multiplet by M^, 
and the mass difference between the inflaton and inflatino by Am^.^ If 
Am^ > m3/2, the decay — + gravitino is kinematically possible. For 



Am^ > few X 7713/2, the partial decay width will be |26, 



Helicity ± - : T 



2 ■ '^^'^+^ 487r Mlmli^Ml 



3 ^ 1 K-^|)' 



'^'^'''^ ^ 2 ■ " 48^ M3Am^M| " ^^'^^ 



We can estimate the abundance of produced gravitinos with the help of total 
inflaton decay rate = ^/{gJ^I^yn^^J9^T^/ My> and the dilution factor due 
to flnal entropy release which is given by 3Tr/4M<^. 

^FuU thermal equilibrium is indeed achieved very rapidly, provided that inflaton decay products 
have interactions of moderate strength. For details on thermalization, see |4^ . 

■^The mass difference between the inflaton and inflatino Am^ is in general different from that 
between the standard model fermions and sfermions m. As an example, consider the case where 
the soft breaking (mass)^ of both the inflaton and sfermion fields is fh? . We will then have Am^ ~ 
{rh? /2M^ <^ fh \i fh <^ M^, while Am^ ~ m when m > M^. 
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If Am^ <^ M^, inflaton decay gives rise to 
Heliaty±i:m ^ / A^VV A™ 

2 V S / reh V "^3/2 / 



Helicity ± - 



3 /7^3/2\ / ^fTT-l 



2 V S /reh VTrMp 



TrMp 
228.75 



228.75 



n 1/2 



1.6 X 10-2 



1/2 

1.6 X 10"^ (2.^ 



An interesting point is that drops out of the calculation, and hence the final 
results in Eq. ( |2.8| ) have no explicit dependence on M^. If < Am^, the 
gravitino abundance will be smaller than that in ( p.8|) by a factor of 16. Note 
again that for Am^ > m3/2 gravitinos of both helicities have approximately the 
same abundance. 

Since from (|2.8| ) ^3/2/5 is inversely proportional to Tr, gravitino production 
from inflaton decay becomes more efficient at lower reheating temperature. 
The reason is that a smaller Tr means a smaller total decay rate T^, while the 
partial decay width ( |2.7| ) is independent from Tr. Therefore decreasing Tr, 
while suppresses the production from thermal scatterings ( p.4|) and sfermion 
decays (|2.6| ), can actually enhance the overall production of gravitinos. Ob- 
viously gravitino production from inflaton decay reaches saturation when the 
partial decay width equals the total decay rate F^. In this case all infiatons 
decay to inflatino-gravitino pairs, and the subsequent decay of inflatinos will re- 
heat the Universe. In consequence, one gravitino will be produced per inflaton 
quanta, resulting in a gravitino abundance (^3/2 = 3Tr/4M0. 

Gravitino production in two-body decays of the inflaton will be kinematically 
forbidden if m3/2 > Am^. However, the inflaton decay inevitably results in 
gravitino production at higher orders of perturbation theory, provided that 
> 1713/2 The leading order contributions come from the diagrams de- 
scribing the dominant mode of inflaton decay with gravitino emission from the 
inflaton, its decay products and the decay vertex. The partial width for infla- 
ton decay to gravitinos is in this case ~ (M^/Mp)'^ which, after taking into 
account of the dilution factor, leads to [713/2/ s) ~ (TjiM^/Mp). If we impose 
the bound on the inflaton mass < 10^^ GeV from the CMB for a simple 
chaotic type inflation model, and if ^3/2 > ^Tn^, gravitino production from 
inflaton decay is subdominant compared to that of thermal scatterings 
and hence can be neglected.'' 



*A similar process, non-thermal production of gravitons from inflaton decay, can become impor- 
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Non-perturbative production: 

We note that besides various perturbative production mechanisms, both of 
the hehcity states can be excited non-perturbatively during the coherent os- 
cillations of the inflaton. This was first discussed in and then elabo- 



rated in |30|. Right after inflation the helicity ±1/2 component, i.e. the 
Goldstino, is essentially the inflatino. For simple models with a single chi- 
ral superfield, it was shown that this component can be produced abundantly 
("-3/2/s) < (Tr/M<^) The reason is that its couplings, given in Eqs. p.l| ) 
and (|2.3|) , are not necessarily Mp-suppressed (contrary to the helicity ±3/2 
states). However, as explicitly shown in |^, it also decays quickly along with 
the inflaton through derivative interactions, and hence poses no danger. Realis- 
tic models include at least two chiral superfields such that the inflation sector is 
different from the sector responsible for SUSY breaking in the vacuum. In these 
models also most of the spin- 1/2 fermions produced during inflaton oscillations 
decay in form of inflatinos, provided that the scales of inflation and present day 
SUSY breaking are sufficiently separated [^]. The helicity ±3/2 components of 
the gravitino have Mp suppressed coupling all the time. In consequence, they 
are produced less abundantly (11.3/2/3) < (M^/Mp) (T^/Mp) pOf, compared 
to the direct decay of the inflaton, thermal scatterings and sfermion decays. 
We will therefore ignore the contribution from non-perturbative production of 
gravitinos in the following. 

To summarize, the total gravitino abundance is given by 

\ S / V S /sea V S /dec V S / reh 

As long as 7713/2 < Tr, gravitinos are always produced in thermal scatterings of 
gauge and gaugino quanta. In addition, sfermion and inflaton decays also contribute 
to gravitino production if m3/2 < < 7r and m3/2 < Am^ respectively. Then it 
turns out from (^), ( ^.6[ ) and (^) that sfermion decays will be the dominant source 



of gravitino production unless Tr, > 1.2 x (^m^ / m^/^^ and/or Am^ > 0.37m. Hence, 
for m3/2 < m < Tr, the most important contribution in (|2.9| ) in general comes from 
sfermion decays (see also footnote 3 on page 6). 



tant in models with extra dimensions [Q . The large multiplicity of the Kaluza-Klein modes of the 
graviton can in this case easily overcome the suppression factor (M^/Mp)^. 
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Gravitino Decay 



Stable gravitino: 

First, we briefly consider the case for stable gravitinos. If the gravitino is the 
LSP, and -R-parity is conserved, it will be absolutely stable. Its total abundance 
(including both helicity ±1/2 and ±3/2 states) will in this case be constrained 
by the dark matter limit ^3/2/?-^ < 0.129, leading to 

_ /I npV\ 

(3.1) 



This implies that the individual contributions from Eqs. ( |2.4| ), ( |2.tj| ) and 
should respect this bound. As an example, consider the case with 777,3/2 = 10 
GeV and Mg ^ 1 TeV. This results in the constraints Tr < 5.5 x 10*^ GeV, 
m < 33 TeV and Am^ < 140 TeV. 

Unstable gravitino: 

An unstable gravitino decays to particle-sparticle pairs through the couplings 
in ( |2.1|) , and the decay rate is given by p9| 

(AT \ JJl'^ 

where Ng and Nf are the number of available decay channels into gauge-gaugino 
and fermion-sfermion pairs respectively. The gravitino decay is completed when 
H ~ r3/2, when the temperature of the Universe is given by 



3/2 



10.75 



.105 GeV/ 

Here g^{T^j2) is the number of relativistic degrees of freedom at T3/2. If 7773/2 < 
50 TeV, gravitinos decay during or after BBN and can ruin its successful 
predictions for the primordial abundance of light elements . If the gravitinos 
decay radiatively, the most stringent bound (773/2/5) < 10~^^ — 10~^^ arises for 
7773/2 — 100 GeV — 1 TeV |^. On the other hand, much stronger bounds are 
derived if the gravitinos mainly decay through hadronic modes. In particular, 
a branching ratio ~ 1 requires that (713/2/5) < lO-^'' - 10-^5 in the same 
gravitino mass range 

To give a numerical example, consider the case when 7/23/2 — 1 TeV. The 
abundance of a radiatively decaying gravitino is in this case constrained to be 



(^3/2/s) < 10-12 |4|_ Thgn Eqs^ (^^^ and Q result in the bounds 

Tr < 10^° GeV, m < 203 TeV and Am^ < 1.1 x 10^ GeV, respectively. If a 
TeV gravitino mainly decays into gluon-gluino pairs, which will be the case if 
m3/2 > Mg, we must have (^3/2/5) < lO^"*^^ This leads to much tighter 
bounds Tr < 10^ GeV, m < 9.4 TeV and Am^ < 11 TeV. 

3.1 Decay of Supermassive Gravitinos and Dark Matter Abundance 

We now turn to supermassive gravitinos with a mass m3/2 > 50 TeV. If one insists on 
a successful supersymmetric gauge coupling unification, the gaugino masses should 
be below 100 TeV. This implies that the gravitino will not be the LSP. The decay 
of supermassive gravitinos happens sufficiently early in order not to affect the BBN. 
Nevertheless, their abundance can still be constrained due to different considerations. 
Gravitino decay produces one LSP per gravitino. This non-thermal component may 
exceed the dark matter limit if the decay happens below the LSP freeze-out temper- 
ature Tf. The freeze-out temperature is given by |jl4 



28 + In 



1 TeV 



10- 



86.25 



1/2" 



(3.4) 



where is the LSP mass and we have parameterized the non-relativistic x annihi- 
lation cross-section as 

{(TxVrel) = —■ (3.5) 



m 



Note that neutralinos reach kinetic equilibrium with the thermal bath, and hence 
become non-relativistic, very quickly at temperatures above MeV [^. The exact 
value of c depends on the nature of x and its interactions. For Bino-like LSP, c can 
be much smaller than for a Wino- or Higgsino-like one. When sfermions are much 
heavier than the neutralinos, c = 3 x 10~^ for a Higgsino LSP and c = 10"^ for a 
Wino LSP (including the effects of co-annihilation) [pl[| . 

Gravitino decay occurs after the LSP freeze-out if T3/2 < T{, which translates 
into an upper bound on the gravitino mass 



m 



■3/2 



< 



1 TeV J 



y/3 






86.25 



1/6 



4.3 X 10^ GeV. 



(3.6) 



This implies that for = 1 TeV, gravitinos with a mass 7713/2 < 4 x 10^ GeV 
decay when thermal annihilation of LSPs is already frozen. This decay produces one 
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LSP per gravitino. The dark matter hmit f^^^^ — 0-129 constrains the total LSP 
abundance to obey 



< 4.6 X 10 



-10 



1 GeV^ 



(3.7) 



The final abundance of LSPs produced from gravitino decay depends on their anni- 
hilation rate. The rate of annihilation of non-relativistic LSPs is given by 



If > r 



3/2, annihilation will be efficient and reduce the LSP abundance to 



(3.8) 



72 

~ 41.58 

s 



86.25 



1/4 



(m3/2Mp)V2- 



(3.9) 



Otherwise, gravitino decay contributes an amount n3/2/s to the LSP abundance. 
Having a large abundance of gravitinos, i.e. (^3/2/^) > 4.6 x 10^^° (1 GeV/m^), is 
therefore potentially dangerous and requires special attention. 

The condition for efficient annihilation of LSPs whose abundance {n^/ s) > 4.6 x 
10^^° (1 GeV/m^) at the time of gravitino decay translates into a lower bound on 
the gravitino mass 



m3/2 



> 



'10-2" 


2/3 


86.25 


1/6 


c 






U TevJ 



2 X 10^ GeV. 



(3.10) 



If m3/2 is in the window given by ( p.6|) and ( |3.10| ), the final abundance of non-thermal 
LSPs will be given by ( p.9| ). It satisfies the dark matter limit ( |3.71 ), and can account 
for the CDM for those values of and ^3/2 which saturate the inequality in (|3.10| ). 
Eqs. ( p.6|) and ( p.lO[ ) can be simultaneously satisfied only if 



' C ' 


1/2 


'fl-* (7^3/2)" 


U0-2J 




86.25 



1.8 TeV. 



(3.11) 



As a matter of fact, this is also the condition such that thermal abundance of LSPs 
at freeze-out respects the dark matter limit. It is not surprising as T3/2 = Tf when 



the inequality in ( |3.11| ) is saturated. For the saturation value of thermal LSP 
abundance gives rise to ^^^"^ ~ 0.129. For smaller the thermal component is not 
sufficient, while for larger thermal LSPs overdose the Universe. 

For values of m-^ respecting the bound in (|3.11|) , there always exists a window 
for m3/2 such that gravitinos decay after the freeze-out while at the same time non- 
thermal LSPs efficiently annihilate and their abundance respects the dark matter 
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limit ( p.7|) . In this mass window the abundance of thermal LSPs is too low to 
account for dark matter. On the other hand, non-thermal dark matter will be a viable 
scenario when the inequality in Eq. ( p.lO| ) is saturated (for non-thermal production 
of LSP dark matter from gravitino decay, see also |]2T], 0). The gravitino mass 
window becomes narrower as increases. It will eventually disappear when 
reaches the upper bound in (|3.11| ). For the canonical choice of c = 10"^, the gravitino 
mass window shrinks to a single point 77^3/2 = 6.3 x 10^ GeV at the saturation value 
= 1.8 TeV. For larger LSP masses it is necessary that gravitinos which decay after 
the freeze-out are not overproduced, i.e. that (^3/2/5) < 4.6 x 10^^° (1 GeV/m^). 
Otherwise, gravitinos must decay above the freeze-out temperature, i.e. the opposite 
inequality as in ( p. 61 ) must be satisfied. Then gravitino decay does not affect the 
final LSP abundance as T3/2 > Tf . However, for masses violating the bound in ( p.llj ) 
thermal LSPs overdose the Universe. Therefore a viable scenario of LSP dark matter 
in this case requires late entropy generation. 

3.2 Gravitino Non-domination 

We now consider the constraints from dark matter abundance on gravitino decay 
in more detail. Assuming that there is no other stage of entropy generation, the 
Universe will remain in the radiation-dominated phase after reheating. During this 
period the scale factor the Universe increases as a oc H~^^^. Gravitinos become 
non-relativistic at ^ 

H,,, ^ (^^^ H^, (3.12) 

where Hp denotes the expansion rate when (most of the) gravitinos are produced and 
Ep is the energy of gravitinos upon their production. If thermal scatterings are the 
main source of gravitino production. Hp ~ T^/Mp and Ep ~ Tr. On the other hand, 
if sfermion decays dominate gravitino production. Hp ~ rn? /M-p and Ep = m/2. 
Finally, if most of the gravitinos are produced in infiaton decay. Hp ~ T^/Mp and 
Ep ^ Am^. 

For H < Hp the energy density of the gravitinos is redshifted oc a^'^, compared 
to oc a^"^ for radiation. Initially the gravitino energy density is P3/2 = "^3/2 -^'P? while 
the energy density in radiation is prad = (7r^/30) (7*(Tp)Tp . Gravitinos will dominate 
when p3/2/prad is compensated by the slower redshift of P3/2. This happens at 



16 ^^3/2^2 fE, 

'dom 



9 \ s J \Tp 



rr _oq ^?*(T,^1^/' /-.ox2m2 

I -'J non O.C/ 



228.75 



n3/2Y "'3/2 . . 
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where we have used s = (27r^/45) (7*(Tp)Tp. Here g*(Tp) is the number of relativistic 
degrees of freedom in the thermal bath at the temperature Tp when gravitinos are 
produced. Gravitino non-domination therefore requires that Hdom < ^3/2, i-e. that 
gravitinos decay while their energy density is subdominant. This translates into the 
bound 



^3/2 ^ 



228.75 



1/4 



_ "^3/2 ^ 

on the gravitino abundance. 

It is seen that for 7/13/2 > 50 TeV, gravitinos will dominate if {n3/2/s) > 1.7 x 
10~^. Thermal scatterings alone, see ( p.4| ), can yield such large abundances for 
extremely large reheating temperatures Tr > 10^^ GeV (we consider Mg < 100 TeV 
here). Therefore they do not lead to gravitino domination in general. The sfermion 
and inflaton decays, however, can produce a sufficiently large number of gravitinos 
for much lower Tr. As mentioned earlier, see the discussion after Eq. (|2.9| ), sfermion 
decays are usually the dominant source of gravitino production when Tr > fh. We 
therefore concentrate on the sfermions here. 

Sfermion decays, see (|2.6| ), will not lead to gravitino domination if 



Here we have taken g^.{m) = 228.75 and = 46 in (2.6). 

A successful scenario with gravitino non-domination should take the constraints 
from LSP production into account. Fig. (1) depicts different regions in the m — 7713/2 
plane for the choice c = 10~^ and = 100 GeV. Above the solid line gravitinos 
dominate, i.e. the opposite inequality as in ( |3.15|) is satisfied, and hence excluded. 
The region between the solid and dashed lines is defined by 



105 GeV J V105 GeV J 



3.4 X 10*^ GeV < m < . Z 1.3 x 10« GeV . 



IGeVV^V ^3/2 A2/%,..,,7n„...^ . ^ "^3/2 ^ 

(3.16) 

In this region {n^/2/s) > 4.6 x 10^^° (1 GeV/?77^), thus the density of LSPs produced 
in gravitino decay exceeds the dark matter limit. The dotted and dot-dashed vertical 
lines correspond to Eqs. ( ^.6] ) and ( |3.10D respectively. The regions in black color are 
excluded since either gravitinos dominate or gravitino decay produces too many LSPs 
which do not sufficiently annihilate. In region 1 gravitinos decay after the freeze- 
out but efficient annihilation reduces the abundance of produced LSPs below the 
dark matter limit. Gravitino decay occurs before the freeze-out in region 2 and 
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does not affect the final LSP abundance. Below the dashed line sfermion decays 
do not overproduce gravitinos. In fact, below the m = m3/2 line such decays are 
kinematically forbidden altogether. In this part of the m — m^j2 plane one has 
to worry about thermal scatterings though. If Tr > 4.6 x 10^° (1 GeV/m^) GeV, 
scatterings will overproduce gravitinos. Therefore regions 3, 4 and 5 will not be 
acceptable in this case, due to inefficient LSP annihilation. 

For the values of c and chosen in this plot, thermal abundance of LSPs at 
freeze-out is too small to account for dark matter. Non-thermal LSP dark matter 
from gravitino decay will in this case be a viable scenario along the dot-dashed line. 




Figure 1: Constraints from LSP production in gravitino non-domination case for c = 10~^ 
and = 100 GeV. The solid and dashed lines represent the upper and lower limits in 
Eq. (3.16) respectively. The dotted and dot-dashed lines correspond to Eqs. ( p. 61 ) and ( 3.1C| ) 
respectively. The regions in black color are excluded. The solid red line is given by 
m = m3/2. 



3.3 Gravitino-dominated Universe 

Gravitinos eventually dominate the energy density of the Universe if T^j2 < -f^dom? 
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Figure 2: Parameter constraints from LSP production in gravitino-dominated case for 
c = 10^^ and = 100 GeV. The dotted and dot-dashed hnes correspond to Eqs. 
and ( |3.10| ) respectively. The region in black color is excluded. 

which happens for 



"-3/2 



> 



228.75 



1/4 



m3/2 \ 1/2 



2.4 X 10" 



(3.17) 



MTp)\ V 105 GeV y 

The scale factor of the Universe a oc H~^^^ in the interval 1^ H < ifdom- Grav- 
itino decay will then increase the entropy density by the factor d = {g^:{Ta,)T^ / g^:(T\j)T^) 
Here Ta, Tb denote the temperature of the thermal bath before and after gravitino 
decay, respectively, while g^iT^^), g^iT^,) are the number of relativistic degrees of free- 

dom at Ta and Tb, respectively. Note that d = ig^,{Ta)pl^2/ 9*i'^b)PR) ? with P3/2 
and pr being the energy density in the gravitinos and radiation, respectively, at the 
time of gravitino decay. The dilution factor d is therefore given by 



d 



H. 



dom 



3/2 



1/2 



1/4 







(^3/2/s) 


_228.75_ 


V "^3/2 / 


2.4 X 10-8 



(3.18) 



^To be more precise, the dilution factor is given by 1 + d. The two definitions are essentially 
equivalent when d <C 1. Obviously there is no gravitino dominaiton, and hence no dilution, when 
d < 1. 
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Here we have taken ((7*(Ta)/(7*(Tb)^^'^ ~ 1, which is a good approximation since in a 
wide range 1 Mev < T < m the number of relativistic degrees of freedom g^, changes 
between 10.75 and 228.75. To be more precise, the dilution factor is given by 1 + c?. 
Obviously for (i < 1 there is no gravitino domination, and hence no dilution. As 
mentioned before, sfermions are usually the main source of gravitino production for 
> fn, and hence we concentrate on them here. 
Gravitinos produced in sfermion decays dominate the Universe if 

(l^G^^)'^' ^ GeV < m < Tk , (3.19) 
in which case the dilution factor is given by 

V "3/2 J 1, 1.3 X lO" GeV ) ■ 

Gravitino decay dilutes the existing LSP abundance by a factor of d, while at the 
same time producing one LSP per gravitino. Note that (^3/2/5) < 10^^ as gravitinos 
at most reach thermal equilibrium when m ^ Eq- ( |3.18|) then implies that 

d < 5.9 X 10^. It is evident from ( |3.17| ) that the abundance of non-thermal LSPs 
upon their production from the decay of supermassive gravitinos exceeds the dark 
matter limit by several orders of magnitude. Therefore any acceptable scenario 
with gravitino domination requires that the LSP annihilation be efficient at the 
temperature T3/2, i.e. that ( |3.10| ) be satisfied. 

Fig. (2) depicts different regions in the m— 7713/2 plane for the choice c = 10^^ and 
= 100 GeV. Gravitinos dominate in the region above the solid line corresponding 
to Eq. ( |3.19|) . Therefore regions below this line are irrelevant. To the left of the 
dotted line, which represents Eq. (|3.6|) , gravitino decay occurs after the freeze-out. 
To the right of the dot-dashed line, corresponding to Eq. ( |3.10| ), LSPs produced from 
gravitino decay efficiently annihilate. The region in black color is excluded due to 
inefficient LSP annihilation. Gravitinos decay after the freeze-out in region 1, but 
the final abundance of LSPs respects the dark matter limit. On the other hand, 
T3/2 > Tf in region 2, and hence gravitino decay does not affect the LSP abundance. 
For the values of c and chosen here, thermal abundance of LSPs at freeze-out is 
too low to account for dark matter. However, non-thermal LSP dark matter from 
gravitino decay is successful along the dot-dashed fine. 

In passing we note that the same discussions also apply when Tr < in. In this 
case, however, the inflaton decay will be the main source of gravitino production 
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as sfermions are not excited by the thermal bath. Constraints from efficient LSP 
annihilation then lead to plots similar to those in Figs. (1), (2), with the m — 7713/2 
plane replaced by the Am^ — m3/2 plane. 

Entropy release from gravitino decay also dilutes any previously generated baryon 
asymmetry. This implies that baryogenesis should either take place after gravitino 
decay, or generates an asymmetry in excess of the observed value by the dilution 
factor given in ( p.l8| ). It is seen from ( p.3| ) that T3/2 < 100 GeV unless m3/2 > 10^ 
GeV. This implies that successful baryogenesis after gravitino decay will be possi- 
ble only if gravitinos are extremely heavy. Therefore the more likely scenario in 
a gravitino-dominated Universe is generating a sufficiently large baryon asymmetry 
at early times which will be subsequently diluted by gravitino decay. We will dis- 
cuss leptogenesis and baryogenesis, and the effect of gravitinos in detail in the next 
Sections. 

One might invoke an intermediate stage of entropy release by the late decay 
of some scalar condensate (beside inflaton) to prevent gravitino domination. We 
shall notice, however, that any such decay will itself produce gravitinos with an 
abundance which is inversely proportional to the new (and lower) reheating temper- 
ature, see (|2.8|) . ^ This implies that any stage of reheating, while diluting gravitinos 
which are produced during the previous stage(s), can indeed produce more graviti- 
nos. Therefore entropy generation via scalar field decay is in general not a helpful 
way to avoid a gravitino-dominated Universe. 

One comment is in order before closing this subsection. In both of the gravitino 
non-domination and domination scenarios, having an LSP abundance in agreement 
with the dark matter limit constrains its mass through Eq. (|3.11 ). Heavier LSPs 
overdose the Universe in one way or another. If T3/2 > Tf, gravitino decay will be 
irrelevant but thermal abundance of LSPs will be too high. If T3/2 < Tf, gravitino 
decay can in addition make an unacceptably large non-thermal contribution. In case 
of gravitino domination gravitino decay dilutes thermal LSPs. However, according 
to (|3.17|) , the decay itself overproduces non-thermal LSPs which will not sufficiently 
annihilate. Therefore gravitino domination cannot rescue a scenario with thermally 
overproduced LSPs. Indeed, for ^ 1 gravitinos should never dominate the Uni- 
verse. The problem can be solved if Tr < Tf, or if another stage of entropy release 
below Tf dilutes thermal LSPs. In both case, however, reheating can overproduce 



^If the scalar field does not dominate the Universe, the expression in (2.8) should be multiplied 
by the fraction r of the total energy density r which it carries. 
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gravitinos and the subsequent gravitino decay may lead to non-thermal overproduc- 
tion of LSPs. If -R-parity is broken, the LSP will be unstable and its abundance will 
not be subject to the dark matter bound. Obviously its thermal and/or non-thermal 
overproduction will not pose a danger in this case. 

3.4 Solving the Boltzmann Equation 

The entropy generated by gravitino decay can be estimated from Eq. (|3.20|) . However, 
the evolution of gravitinos and relativistic particles can be followed directly by solving 
the Boltzmann equation. Assuming that gravitinos are non-relativistic at decay and 
that g^: is constant during the decay process, the Boltzmann equations for gravitinos 
take the form 

P3/2 = -r3/2P3/2 - 3iJp3/2 (3.21) 

and 

f=-HT + Ty2^^T^. (3.22) 
These should be solved together with the Friedmann equation 



3M2 



7r2 



P3/2 + g*iT)^T'^ 



(3.23) 



In Fig. (^ we show several examples of solving the Boltzmann equation for different 
initial conditions, always assuming the = 10.75 during decay. The dilution factor 
d can then be found from 

where a[ and Ti denote the initial values of the scale factor and the radiation tem- 
perature respectively. The bottom panel of Fig. (H) shows the dilution factor, and 



when gravitinos dominate it agrees quite well with Eq. (3.20). 



4. Leptogenesis 

The baryon asymmetry of the Universe (BAU) parameterized as rj^ = [ub — n^)/ s, 
with s being the entropy density, is determined to be 0.9 x 10^^° by recent analysis of 
WMAP data p|. This number is also in good agreement with an independent deter- 
mination from primordial abundances produced during BBN Three conditions 
are required for generating a baryon asymmetry: B— and/or L— violation, C— and 
CP— violation, and departure from thermal equilibrium [Q. Since B + L- violating 
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Figure 3: Plot of temperature, density and entropy increase as functions of time for 
"^3/2 = 100 TeV. Curves are for in^/2/s) = 10~^^, 10^^, 10"^, 10^'', and 10^^ in increasing 
order. 



sphalerons transitions are active at temperatures 100 GeV <T< 10^^ GeV any 
mechanism for creating a baryon asymmetry at T > 100 GeV must create a, B — L 
asymmetry. The final asymmetry is then given hy B = a{B — L), where a = 28/79 
in case of SM and a = 8/23 for MSSM [g. 

Leptogenesis postulates the existence of RH neutrinos, which are SM singlets, 
with a lepton number violating Majorana mass M^. It can be naturally embedded 



in models which explain the light neutrino masses via the see-saw mechanism [51 



A lepton asymmetry can then be generated from the out-of-equilibrium decay of the 
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RH neutrinos into Higgs bosons and light leptons, provided CP— violating phases 
exist in the neutrino Yukawa couplings 



53 



|. The created lepton asymmetry 
will be converted into a baryonic asymmetry via sphalerons processes. 

The on-shell RH neutrinos whose decay is responsible for the lepton asymmetry 
can be produced thermally via their Yukawa interactions with the standard model 
fields and their superpartners [^], for which Tr > Mi ~ 10^ GeV, |^ |56|, |5^, 
or non-thermally for which Tr < M^r, see ^ Non-thermal leptogenesis can 
also be achieved without exciting on-shell RH neutrinos |6^. In supersymmetric 
models one in addition has the RH sneutrinos which serve as an additional source for 



leptogenesis . The sneutrinos are produced along with neutrinos in a thermal bath 
or during reheating, and with much higher abundances in preheating . There are 
also additional possibilities for leptogenesis from the RH sneutrinos some of which 



rely on soft SUSY breaking effects |6a, pE 37, 38 



4.1 Thermal Leptogenesis 

Let us concentrate on the supersymmetric standard model augmented with three 
RH neutrino multiplets in order to accommodate neutrino masses via the see-saw 
mechanism The relevant part of the superpotential is 



(4.1) 



where N, H^, and L are multiplets containing the RH neutrinos and sneutrinos 
A^, the Higgs field giving mass to, e.g., the top quark and its superpartner, and 
the left-handed (s) lepton doublets, respectively. Here hij are the neutrino Yukawa 
couplings and we work in the basis in which the Majorana mass matrix is diagonal. 
The decay of a RH (s) neutrino with mass Mi (we choose Mi < M2 < M3) results in 
a lepton asymmetry per (s) neutrino quanta e^, given by 



_L_l_5:,.(,hhtyV(|) 



(4.2) 



with 54 
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X 



1 + x 



X 



(4.3) 



The first and second terms on the right-hand side of Eq. (^4.3| ) correspond to the one- 
loop self-energy and vertex corrections, respectively. Assuming strongly hierarchical 
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RH (s)neutrinos and an 0{1) CP— violating phase in the Yukawa couphngs, it can 



be shown that 71 



3 Mi(m3-mi) 

where mi < m2 < 77^3 are the masses of hght, mostly left-handed (LH) neutrinos. 
For a hierarchical spectrum of light neutrino masses (mi <C m2 <C rn^), we then have 

To obtain this, we have used m3 — mi ~ m3 ~ 0.05 eV (as suggested by atmospheric 
neutrino oscillation data) and (Hy) — 170 GeV. 

If the asymmetry is mainly produced from the decay of the lightest RH states, 
after taking the conversion by sphalerons into account, we arrive at 

where we have assumed maximal CP- violation. Here k is the efficiency factor ac- 
counting for the decay, inverse decay and scattering processes involving the RH 
states m, 13]- 

The decay parameter K is defined as 

where 

i 

is the decay width of A^i and Ni. One can also define the effective neutrino mass 

Ml ' (^-^^ 

i 

which determines the strength of Ni and A^i interactions, with the model-independent 
bound mi < mi 1173 . 



case(l): 

If < 1, corresponding to mi < 10~^ eV, the decay of RH states will be out-of- 
equilibrium at all times. In this case the abundance of RH states produced via 
Yukawa interactions does not reach the thermal equilibrium value. The lepton 
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number violating scatterings can be safely neglected. Hence this case is called 
the weak washout regime. The efficiency factor is k ~ 0.1 when rhi = 10~^ eV. 
Generating sufficient asymmetry then puts an absolute lower bound few x 10^ 
GeV on Mi, and Tr, > Mi will be required in this case [BS . 



• case (2): 

In the opposite limit K > 1, Ni and A^^i will be in thermal equilibrium at 
temperatures T > Mi. In particular, the efficiency of inverse decays erases any 
pre-existing asymmetry (generated, for example, from the decay of heavier RH 
states). This regime is called of strong washout. We note that this regime in- 
cludes the entire favored neutrino mass range mgoi ^ fhi < matm- The inverse 
decays keep the RH (s)neutrinos in equilibrium for sometime after T drops 
below Ml. The number density of quanta which undergo out-of-equilibrium 
decay is therefore suppressed and reduces the efficiency factor. Successful lep- 
togenesis in the range (mgoi, ""^atm) requires that 10^° GeV < Mi < 10^^ GeV 
while, due to the efficiency of inverse decays, Tr can be smaller than Mi by 
almost one order of magnitude |55|. The efficiency factor n in this window 
varies between few x 10~^ and few x 10~^. 

It is possible to obtain (approximate) analytical expressions for the efficiency 
factor K. In the strong washout regime, the final efficiency factor is maximal when 
AL = 2 scatterings among the LH (s)leptons can be neglected. So long as the 
scatterings can be neglected, the efficiency factor k is independent from the lightest 
RH (s) neutrino mass Mi and K(mi) is given by ^ Q ^ 

..lO-f^V'. (4.10) 



mi 

Scatterings cannot be neglected for very large Mi or large light neutrino masses. In 



this case the efficiency factor can be approximated as ||55| : 



K{mi,Mi,m^) = K{mi)e -b'^io"' gcvAi cv; ^ (4.11) 

where u ~ 0.186, zb = Mi/T-q ~ few, with Tb the temperature at which most of the 
lepton asymmetry is produced, is the sum over the squares of the light neutrino 
masses. The efficiency factor /t(mi) is given in Eq. ( [4.10|) . For hierarchical light 



neutrinos, rn? = m\ ~ Am^^„ ^ 2.2 x 10 ^ eV^, where Am^j^ is the mass squared 
difference which controls the oscillations of atmospheric neutrinos. Washout effects 
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due to scatterings then become important for Mi ~ 10^^ GeV. For quasi-degenerate 
neutrinos, nii ~ m2 — ma, larger values of fh are possible and n will be exponentially 
suppressed even for much smaller values of Mi. 

For negligible AL = 2 scatterings, the baryon asymmetry is, in the case of 
maximal decay asymmetry, given by 

. 0.9 X 10-'« f f ^) " ( % ) . (4.12) 

The baryon asymmetry is proportional to Mi. The usual bound on the reheating 
temperature Tr < 10^° GeV, given for m^i2 < 50 TeV, does not apply to super- 
massive gravitinos. Hence, RH (s)neutrinos with masses Mi ^ 10^'' GeV can be 
thermally produced in the early Universe. This can lead to larger amounts of baryon 
asymmetry generated by thermal leptogenesis, with respect to the standard sce- 
nario. The lowest reheating temperatures required for thermally producing these 
RH (s)neutrinos is ~ Mi/ few <C 10^^ GeV. Note that the initial assumption 
of gravitino non-domination is satisfied if gravitinos are dominantly produced by 
thermal scatterings. 

4.2 Effects of the Gravitino on Thermal Leptogenesis 

Thermal leptogenesis completes when T ~ Mi/ few ||5^, |5^. Eq. ( p.3|) then implies 
that gravitino decay takes place after leptogenesis unless they are extremely heavy: 

On the other hand, for 7713/2 > 50 TeV, gravitino decay occurs at a temperature 
^3/2 > 6.8 MeV which is compatible with a successful BBN. Therefore both sce- 
narios of gravitino domination and non- domination are in agreement with the BBN 
constraints. Nevertheless the effect of gravitino decay on the final baryon asymmetry 
need to be taken into account. We consider both scenarios of gravitino domination 
and non- domination. 

• Gravitino non- domination: 

The condition for gravitino non-domination is given in Eq. ( glil) . There will 



be practically no dilution by gravitino decay in this case, and thermal lep- 
togenesis should generate rj^ — 10^^° according to Eq. ( [4.6| ). Obtaining suf- 
ficient asymmetry in both of the weak and strong washout regimes requires 
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that Ml > Tr > 10^ GeV, and leptogenesis completes when T ~ Mi > 10^ 
GeV [p5| , p7[] . Sfermions with a mass m < 10^ GeV certainly reach thermal 
equilibrium and, for m^i2 < m, their decay will produce gravitinos according 
to Eq. 



If m3/2 > gravitino production from sfermion decays is kinematically for- 
bidden. Scatterings of gauge and gaugino quanta in the thermal bath will 
nevertheless produce gravitinos so long as 1713/2 ^ ^r- Late time domination 
of gravitinos thus produced requires that the reheating temperature Tr > 10^"^ 
GeV, see Eq. ( p.4|) . However, gravitinos can be overproduced for much smaller 
Tr. For = 100 GeV (1 TeV), gravitino decay results in LSP overproduc- 
tion when Tr > 3 X 10^° (3 x 10^) GeV. This indeed occurs for the bulk of 
the parameter space compatible with thermal leptogenesis, particularly in the 
favored neutrino mass window mgoiar < ^hi < matm The condition for 



sufficient annihilation of non-thermal LSPs in this case sets a lower bound on 
the gravitino mass through Eq. ( |3.10| ), independently of whether m3/2 < m 
or 1713/2 > m. Constraints from LSP annihilation, which determine acceptable 
regions of the m — 7/13/2 plane, are summarized in Fig. (1) and the related 
discussion. 

In the case of gravitino non-domination, the produced lepton asymmetry is not 
subsequently diluted by gravitino decays, even if taking place after leptogen- 
esis. As the bound on the reheating temperature Tr < 10^° GeV does not 
apply for supermassive gravitinos, RH neutrinos with masses Mi ^ 10^° GeV 
can be fully thermalized for a sufficiently large reheating temperature. From 
Eqs. ( [4.6| ), ([4.10|) and ( [4.11| ), it follows that the baryon asymmetry is pro- 
portional to the lightest RH neutrino mass. Mi, up to Mi ~ 10^^ GeV, for 
~ Am^^^. For larger values of Mi, the lepton asymmetry is washed out 
by AL = 2 scatterings. Then, the maximum baryon asymmetry is produced 
for Ml ~ 5 X 10^^ zb/uj GeV. Here we have again taken ~ Am^^^. From 
Eq. ( ^.12| ), we notice that for large values of Mi the generated rj^ can be much 



larger than the one required to explain the observations, if the decay asymme- 
try is maximal. Models which implement the see-saw mechanism of neutrino 
mass generation typically assume the conservation of flavor symmetries and/or 
special forms of the Yukawa couplings, in order to explain the low energy neu- 
trino masses and mixing. In many of these models, the decay asymmetry is 
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constrained to be non maximal and Mi larger than the typical values 10^- 
10^° GeV are needed to generate a sufficient baryon asymmetry. Models with 
supermassive gravitinos allow to have reheating temperatures high enough to 
thermally produce such heavy RH neutrinos. In each specific model, a detailed 
analysis is required for establishing the feasibility of successfuU leptogenesis 
and, at the same time, the possibility to explain the low energy neutrino mass 
matrix (for a discussion of CP- violation in specific see-saw models and lepto- 
genesis, see, e.g., Ref. [[T^] ). 

Gravitino- dominated Universe: 

If gravitinos are produced very abundantly, see Eq. (|3.171 ), the Universe will 
become gravitino-dominated. Sfermion decays (which, as mentioned earlier, 
usually dominate over thermal scatterings and infiaton decay) produce such 
abundances of gravitinos when Eq. ( ^.19| ) is satisfied. 

Gravitino decay reheats the Universe to a temperature T3/2, see Eq. (^.31), and 
increases the entropy density by a factor d given in Eq. (p.l8|) . Successful ther- 
mal leptogenesis after gravitino decay will be only possible if gravitinos are 
extremely heavy 7713/2 > 10^^ GeV, see Eq. ( [4.13| ). According to Eq. ( |3.17| ), 
gravitino domination in this case requires that (^3/2/5) > 10~^. It follows from 
Eq. ( |3.19D that sfermion decays yield this only if rh (and Tr) is > 10^^ GeV. In 
addition, non-perturbative production of gravitinos with the necessary abun- 
dance is also questionable. 

Therefore, it is more realistic to consider the opposite situation where lepto- 
genesis occurs before gravitino decay. In this case, due to the entropy release 
by gravitino decay, the generated asymmetry must exceed the observed value 
of ?7b — 10~^° by a factor of d. Eqs. (|3.2CI|) and (|4.6|) then imply that the final 
asymmetry is 

' V 105 GeV/ V ^ / \10^ GeV J ^ ' 

As a specific example, consider leptogenesis in the favored neutrino mass win- 
dow TTisoi < rhi < matm- In this interval, which entirely lies within the strong 
washout regime, the efficiency factor k ~ 10~^, and the reheating temperature 
follows Tr > 0.1 Ml [^. Therefore successful leptogenesis requires that 

M.^f^)'7..-J^V3xlO-«GeV. (4.15) 



"^3/2 



1.3 X 108 GeV 
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In Fig. 1^ we show the value of Mi needed to produce the correct baryon asym- 
metry of rj-B = 0.9 X 10~^° as a function of the hghtest RH (s)neutrino mass 
m3/2 and m. The plot is produced assuming k = 10~^, but since ?7b oc k it 
can easily be rescaled for other values. Without dilution, successful leptogen- 
esis for k, = 10~^ requires that Mi ^ 3 x 10^*^ GeV, see Eq. ([4.6|) . Therefore 
having contours with Mi > 3 x 10^*^ GeV indicates gravitino domination. As 
mentioned earlier, thermal leptogenesis fails for Mi > 10^^ GeV due to the 
erasure of generated asymmetry by AL = 2 scattering processes. This hap- 
pens in the light colored region, and hence excludes it. The overlap between 
Figs. (1) and (3) combines the constraints from leptogenesis and dark matter 
considerations. 

Thermal leptogenesis cannot generate a baryon asymmetry which exceeds 10^^. 
The maximal value is obtained in the (hypothetical) case when RH (s) neutrinos 
have thermal equilibrium abundance, and the efficiency factor k and asymmetry 
parameter |ei| are both 1, see Eq. (|4.6| ). This implies that successful leptogene- 
sis in a gravitino-dominated Universe would be impossible if the dilution factor 
d was larger than 10^. However, since gravitinos can at most reach thermal 
equilibrium, we always have d < 5.9 x 10^, see the discussion after Eq. (|3.2CI| ). 
This ensures that there will be no such case where thermal leptogenesis is 
absolutely impossible in a gravitino-dominated Universe. 



5. Baryogenesis from Supersymmetric Flat Directions 

There are many gauge invariant combinations of the Higgs, squark and slepton fields 
along which the scalar potential identically vanishes in the limit of exact SUSY. 
Within the MSSM there are nearly 300 fiat directions which are both F— and D— fiat 
and conserve R-parity [0 . Soft terms, as well as non-renormalizable superpotential 
terms, lift the fiat directions when SUSY is broken. 

A homogeneous condensate along can be formed along a fiat direction in the 
infiationary epoch, provided that the fiat direction mass m is smaller than the Hubble 
expansion rate during infiation. The condensate starts oscillating coherently when 
the expansion rate if ~ m. During this epoch the infiationary fiuctuations of the 
condensate can be converted to density perturbations [^. The condensate can also 
help an efficient reheating to the SM degrees of freedom |7B|. In addition, it can 
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log(m3/2/GeV) 



Figure 4: The value of log(Mi) (in GeV) needed to produce rjB = 0.9 x 10^^*^ as a function 
of m.3/2 and m in the gravitino-dominated scenario. We have chosen the efficiency factor 
K = 10~^, which is the typical value in the favored neutrino mass window. The dark 
color region does not have gravitino domination, while in the light color region thermal 
lepotgenesis fails due to washout by lepton number violating scatterings. 



excite vector perturbations to explain the observed large scale magnetic field . 
If the condensate carries a non-zero baryon and/or lepton number, then the flat 
direction dynamics can be responsible for baryogenesis via the Affleck-Dine (AD) 
mechanism |]7^ (for a review, see Ref. [0). In the following section we make a 
general discussion on the viability of AD baryogenesis for arbitrarily heavy gravitinos 
and/or sfermions. 

5.1 Late Baryogenesis via the Affleck— Dine Mechanism 

The scalar potential for a flat direction (p (not to be confused with the infiaton field 



in Section II) is given by 

V = (ft' + c„i/') + (^Ar + /..c.) + ^l^r'"-' . (5.1) 

Here and chH"^ are the soft mass^ from SUSY breaking in the vacuum and SUSY 
breaking by the non-zero energy density of the infiaton respectively. Here ch can 
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have either sign, and its sign is also affected by radiative corrections |7^, 7^]- The 
last term on the right-hand side of (|5.1| ), which lifts the fiat direction, arises from 
a non-renormalizable superpotential term (i.e. n > 4) induced by new physics at 
a high scale M. In general M could be a string scale, below which we can trust 
the effective field theory, or M = Mp. SUSY breaking in the vacuum and by the 
infiaton energy density generate A-terms A and aH, respectively, corresponding to 
this non-renormalizable superpotential term. For minimal Kahler terms and in case 
of gravity mediation m ~ A ~ 1713/2 and < c// ~ 1. Depending on the symmetries 
of the infiaton sector, both a ~ 0{1) and a <^ 1 are possible. 
The equation of motion for the fiat direction is given by 

dV 

+ 3i/0+ — = 0. (5.2) 
The evolution is easiest to analyze by the field parameterization 

0=-^¥.e'^ (5.3) 
where (p, 9 are real fields. Then the scalar potential can be written in the form 
V(^. 0) = 1 (ra' + c,H^) / + ^^^^^ ^» + .^<"-" . (5.4) 



Here 

f{d) = \A\ cos{n9 + 9a + 9x) + \a\H cos{n9 + 0, + 9x) , (5.5) 

with 6*^, 9a, 6'a being the angular directions for A, a, A respectively. The baryon/lepton 
number density is given by 

nT,,i^ = ^^{ct>*4>-H*) =(39if\ (5.6) 

where (3 is the baryon/lepton number carried by the fiat direction. At the minima 
of the of potential 



_2 2"/"~2 M"~^ 



YyM^{-/W±[/W'-4(n-l)(m2 + CHi/^)]^/^} , (5.7) 



and n6'mm = (2p + l)n - 9a-9x if \a\H > \A\, while n9mm = (2p + l)n -9a-9x if 
\a\H < (with p = 0, 1, . . . , n - 1). 

The radial field quickly settles at one of the minima given in ( |5.7|) during 
inflation. If \a\H ^ 1^41, the phase field ip9 (since 9 is dimensionless) has a mass of 
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order \a\H and it ends up in one of the discrete minima nOiam = — da — d\- The 
phase field has a mass <^ H if \a\H <^ \A\, and hence it freezes at a random value. 

After infiation the Hubble rate decreases as the Universe expands, and so does 
(/9min- The if field tracks the instantaneous minimum of the potential, so its evolution 
can be qualitatively understood by looking at the evolution of the minimum. Once 
fn? ~ \ch\H'^i the minimum of the potential changes from v^min to ^9 = in a 
non-adiabatic manner. At this time starts oscillating in the radial direction with 
frequency rh. The motion of the phase field, which is necessary for generating a 
baryon/lepton asymmetry, requires the exertion of a torque. If \a\H ~ a non- 
adiabatic change in the position of the minimum from n6'niin = '^ — da — d\^o n9^^^■^■^^ = 
71 — 0A — 0\ generates a torque. If \a\H <^ \A\, the freezing of the phase field at 
a random value generates the torque and leads to its motion towards the minimum 
nOmin = — Ga — d\ if \a\H ^ \A\. The potential along the angular direction will 
quickly decrease due to the redshift of see Eq. ( ^.4|) , once ip starts its oscillations. 
In consequence, starts freely rotating in the angular direction at which time a net 
baryon/lepton asymmetry is generated. 

Based on Eqs. (|5.2|) and ([5. 6]), the baryon/lepton asymmetry obeys the equation 



dV 

nB,L + SHnB,L = -P-g^ . (5.8) 
This can be integrated to give at late times t ^ H^^l (see |jl3[) 

^2(n - 2) sin5 , ^ , 2 /r 

"^'^"^ 3(n-3)(go3ct)^ '^'^°- ^'-'^ 

where v^osc, -f^osc denote the value of ip, H when the condensate starts oscillating. 
Here 5 is a measure of spontaneous CP- violation in the potential and sin 5 ~ 1. 
The baryon to entropy ratio will then be 

^B,L 3Tr riuL 



4pR ^M^Hl/ ^^-^^^ 

We parameterize the A-term as \A\ = 'yih. For gravity- mediated SUSY breaking 
typically m ~ |A| ~ 1713/2, while in gauge-mediated models \A\ ~ 7713/2 <^ rh. In split 
SUSY A <^ m, see the discussion in Ref. [^. Here we consider 7 to be arbitrary, 
and hence rh and |y4| be unrelated. 

Since H^sc ^ m and (/^^-^ ~ 2(""2)/2jvfn-3~/|^|^ 

nB,L n-2 \A\T^f m y/^""') 

' ' , (5.11) 



6(n-3) m2 \^|A|M] 
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where M = Mp is taken. 

For the lowest dimensional non-renormalizable term n = 4, and A ~ C>{1), 
the generated baryon asymmetry is n-B,h/s ~ 0.1{\A\Tii/mMp). If \A\ ~ fh then 

~ 10^ GeV is adequate to generate baryon asymmetry of order 10^^°. If \A\ <^ 
m, one would require even larger Tr. However the most desirable feature of AD 
baryogenesis lies in its flexibility to generate a desirable baryon asymmetry even for 
very low reheating temperatures. For instance, if m ~ 10^ GeV and n = 6, the 
required asymmetry can be generated for Tr ~ 10^ GeV when \A\ = m. 

At late times, i.e. H fh, contributions from SUSY breaking by the inflaton 
energy density are subdominant to soft terms from SUSY breaking in the vacuum. 
If |y4p < 4(n — l)m^, the potential has only one minimum at ip = 0. However, 
another minimum appears away from the origin when 4(n — l)m^ < see (|5.7| ). 
In the AD scenario the field starts at large values of and hence it gets trapped 
in this secondary minimum in the course of its evolution in the early Universe. 
If 4(n — l)m^ < < n^fn?, the true minimum is still located at the origin. 

Tunneling from the false vacuum could still save the situation in this case. However, 
for \A\ > nih, the true minimum will be at 7^ 0. Since flat directions have non- 
zero charge and color quantum numbers, this will lead to an unacceptable situation 
with charge and color breaking in the vacuum. It is therefore necessary to have 
|Ap < 4(n — l)m^ in order to avoid entrapment in such vacuum states. This is the 
case in gravity-mediated and gauge-mediated models, as well as split SUSY. Note 
that the same discussion applies to the soft terms induced by non-zero energy density 
of the inflaton. However, these terms disappear at late times and will be irrelevant 
in the present vacuum. 



Thermal effects: 

According to the potential given in Eqs. ( |5.lD and ( |5.4| ), the fiat direction con- 
densate starts oscillating when if ~ m. However thermal effects from reheating 



may trigger an earlier oscillation and lead to a larger value of Hose ||83' 131 • '^^^ 
inflaton decay (in the perturbative regime) is a gradual process which starts 
after the end of inflation. Hence, even before the inflaton decay is completed, 
the decay products constitute a thermal bath with instantaneous temperature 
T ~ (//T^Mp)^/^ > Tr [^. The flat direction has gauge and Yukawa cou- 
plings, collectively denoted by y, to other fields. Its VEV gives a mass ~ yip to 
these fields. If y(p < T, these fields are excited in the thermal bath and reach 
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Figure 5: The reheating temperature as a function of the scalar mass, given by Eq. ( 5.11j ), 



for successful AD baryogenesis. The cases n = 4 and n = 6 are represented by solid and 
dashed lines respectively. The black and red lines are plotted for |^| = m and |^| = 1 TeV 
respectively. 

thermal equilibrium. This, in turn, results in a thermal correction ~ yT to the 
mass. If yT exceeds the Hubble parameter at early times, i.e. for H ^ fh, 
the condensate starts early oscillations [^. On the other hand, if yip > T, the 
fields coupled to will be too heavy to be excited. They will decouple from 
the running of gauge coupling(s) at temperature T instead, which induces a 
logarithmic correction to the free energy ~ log (T/ty^). This triggers early 
oscillations of the condensate if T'^/'p > H when H ^ fh Note that ac- 
cording to Eq. (|5.10|) a larger value of iJosc results in a smaller baryon/lepton 



asymmetry. 

Refs. |8y, ^ have studied thermal effects for the conventional case with m ~ 1 
TeV. Thermal corrections of the former type will become less important as m 
increases. The reason is that y^""^ oc H, see (|5.7|) , and hence y^p > H will be 
more difficult to satisfy for larger m. Also, since if oc at early times, (T^/yj) 
increases more slowly that H . Therefore thermal corrections of the latter type 
will also become less important when m ^ 1 TeV. 
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5.2 Effects of Gravitino on Affleck-Dine Baryogenesis 

In Eq. ( |5.11| ) no specific assumption is made about the source wliicli relieats tlie 
Universe. It can be either the infiaton decay, as usually considered, or the decay of 
the fiat direction condensate. In case the infiaton decay reheats the Universe Tr < m 
and m < Tr are both possible. However, Tr > m if the fiat direction is responsible 
for reheating the Universe. The energy density in the condensate oscillations is 
fn^ip^ . The flat direction has gauge and Yukawa couplings to other fields through 
which it induces a mass oc ip for the decay products. The condensate will decay 
no later than the time when p < m. Hence, since energy density in the radiation 
is ~ T^, we will have Tr > m in this case. According to Eq. the yielded 

baryon asymmetry is oc Tr. This implies a larger asymmetry for larger values of 
Tr. Having an unacceptably large baryon asymmetry is indeed typical when the flat 
direction condensate has a very large VEV such that it dominates the energy density 



and, subsequently, reheats the Universe [0. Gravitino domination can in this case 
help to dilute the excess of baryon asymmetry. It is interesting that such a solution, 
invoked from the early days of AD baryogenesis |^, can be naturally realized with 
supermassive gravitinos. 

When Tr > m sfermions reach thermal equilibrium after reheating and their de- 
cay will be the dominant source of gravitino production. The condition for gravitino 
domination and the dilution factor from gravitino decay are then given by Eqs. ( p.l9| ) 
and ( |3.2(J| ) respectively. The final asymmetry generated via the AD mechanism in 



case of gravitino domination will then be, see (|5.11 



ne ^ n-2 / m V'""'^ (\A\\ ^ mg/a \5/2 / 1.3 x Ips GeV V 

The results for the two extreme cases |y4| = m and \A\ = 1 TeV are summarized 
in Fig. (6) when n = 4, 6. The combined constraints from baryogenesis and dark 
matter considerations will be included in the overlap of Figs. (1) and (6). 

If Tr -C rh, sfermion quanta will not be excited in the thermal bath. However, 
infiaton decay can in this case result in efficient production of gravitinos according 
to Eq. ([^.81). One can then repeat the same steps to find an expression for the initial 



asymmetry similar to that in (|5.12|) . Such an expression, and plots similar to those 
in Fig. (6), will however depend on Am^ as well as rh and m3/2. This leads to a more 
complicated and model-dependent situation. Moreover, the scenario with gravitino 
domination will be more constrained when Tr <^ m (specially if \A\ <^ fh). The 
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Figure 6: The scalar mass as a function of gravitino mass, given by Eq. ( 5.12 ), for 
successful AD baryogenesis in a gravitino-dominated Universe. The conventions are the 
same as in Fig. (5). 



5.3 Late Gravitino Production from Q-ball Decay 

So far we have assumed that gravitinos are mainly produced in sfermion decays 
(if > fn), or in inflaton decay (if Tr <^ fh). The flat direction condensate 
consists of zero-mode quanta of the sfermions, and hence its decay too can lead 
to gravitino production. If Tr > m, this will be subdominant to the contribution 
from the decay of thermal sfermions. The reason is that the zero-mode quanta have 
at most an abundance (n/s) which is comparable to that of thermal sfermions. If 
Tr -C m, sfermions will not be excited in the thermal bath. As mentioned before, the 
condensate certainly decays no later than the time when ip < m. Even if ~ Mp 
initially, the condition ip < m is satisfied a.t H > m^/M-p. Not that the scale 
factor of the Universe a oc H^^^^ during reheating, in which phase the Universe 
is dominated by inflaton oscillations. Also, the abundance of zero-mode quanta 
(n/s) < {3Tji/4rh) <^ 1. This implies that the condensate contains a much smaller 
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number of quanta which survive (much) shorter than thermal sfermions. Therefore 
gravitino production from the decay of the flat direction condensate will not be as 
constraining as that in the decay of thermal sfermion. In most case, it can be simply 
neglected. 

The above discussion strictly applies to the oscillations of a homogeneous con- 
densate. However, it usually happens that the flat direction oscillations fragment and 
forms non-topological solitons known as Q-balls p3| . These Q-balls can decay much 
later than a homogeneous condensate. Their late decay could then efficiently produce 
gravitinos, even if Q-balls do not dominate the energy density of the Universe. 

To elucidate, let us consider the case when tree-level sfermion masses at a high 
scale M are given by fh. The potential for the sfermions, after taking into account 
of one-loop corrections, reads 



V = m^\<p\' 



|2\ n 



M2 



(5.13) 



where i^' is a coefficient determined by the renormalization group equations, see [11 



8^ . In order to form Q-balls it is necessary that the potential be flatter than |0p at 
large field values, i.e. that K < 0. Loops which contain gauginos make a negative 
contribution oc —^\i2 to with mi/2 being the gaugino mass. On the other hand, 
loops which contain sfermions contribute oc Then K < Q can be obtained. 



provided that 2mi/2 > m In models of gravity-mediated SUSY breaking m ~ 
mi/2, and hence < is obtained for many flat directions. When the spectrum is 
such that m ^ ^1/2, like in the case of split SUSY, there are no Q-balls as > 0. 

The potential can also be much flatter oc ln|0| at large fleld values. This happens 
in models of gauge- mediated SUSY breaking pSI, and can arise from thermal cor- 



rections . The important point in any case is that the scalar fleld proflle within a 
Q-ball is such that the fleld value is maximum at the center yjo and decreases towards 
the surface. This implies that for gi^Q > m, with g being a typical coupling of the (j) 
fleld to other flelds, the Q-ball decays through its surface as decay inside the Q-ball 
is not energetically allowed [^. Note that for a typical gauge or Yukawa coupling 
g(fo > m ii (fo ^ 'fn. The decay rate of Q-ball which contains a (baryonic/leptonic) 
charge Q is in this case given by |86| 

^ < (5.14) 

dt - 1927r2' ^ ' 

where u c:^ in and A = 4:7tRq is the surface area of the Q-ball. For example, consider 
Q-ball formation for the potential given in ( p.l3| ). A Q-ball with total charge Q then 
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has a decay lifetime [B4] 



-Q^(^Vi5^VjHlxlO-.ec, (5.15) 



0.037 \ fh J 
which corresponds to the decay temperature 

nnQ\l/2 / ~ \ 1/2 / in20 \ 1/2 

0.03 \ / m \ I 10 ^ 



Here we have used r^^ ~ (T^/Mp). The total baryonic/leptonic charge Q of a Q-ball 
is given by the multiplication of baryon/lepton number carried by the flat direction 



and the total number of zero- mode quanta inside the Q-ball. This, after using ( |2.5| ), 
leads to 

Helicity ± i : f . f ^ V f ^ V f « ) f 2.5 X 10- . 



2'\ s /Q-ball V^3/2/ VlTeVy 



s 



Obviously s) < {ub/s), since the decay of each quanta inside the Q-ball 

can at most produce one gravitino. This implies that gravitinos produced from the 
decay of Q-balls will not dominate the Universe if the decay generates a baryon 
asymmetry (ub/s) ~ 10^^°, see Eq. ( p.l7| ). However, the situation will be different 
for larger Q-balls which yield an asymmetry ^ 10^^°. If gravitinos from Q-ball decay 
dominate the Universe, they will dilute the baryon asymmetry. The final asymmetry 
will then have the correct size, see Eqs. (|3.18|) and ( |5.17|) , if 

Assuming that the Q-balls do not dominate the energy density of the Universe, we 
must have (ub/s) < (3Td/4m). The necessary condition for gravitino domination, 
after using Eqs. ( |3.17| ), ( |5.16| ) and Eq. ( |5.17| ), is then obtained to be 

For m ^ 1 TeV and m3/2 "C fh, this lower bound on the Q-ball charge is compatible 
with the value in Eq. (|5.18|) required for successful baryogenesis. Note that the Q- 
ball decay temperature may be smaller than the LSP freeze-out temperature Tf. 
In this case Q-ball decay can be dangerous as three LSP per baryon number will 
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be produced. However, this will not lead to problem so long as the gravitino decay 
temperature T3/2 < Td (which is typically the case) and the condition for efficient 
LSP annihilation in Eq. ( |3.10| ) is satisfied. The Q-balls will dominate the Universe if 
they carry a very large charge. The initial baryon asymmetry released by the Q-balls 
then has a simple expression (ne/s) ~ (Td/m). 

The results in Eqs. (|5J[5|) , f ^) , (113), (ISlSD and (|5l9|) are valid for the 
potential given in Eq. ( ^.13| ) . The same steps (though more involved) can be followed 
to obtain similar results for logarithmic potentials as in the case of gauge-mediated 
models. The remarkable point in all cases is that the Q-ball decay lifetime increases 
with its charge, implying a more efficient production of gravitinos from the Q-ball 
decay. This leads to an attractive solution that the large baryon asymmetry released 
by the Q-ball decay can be naturally diluted by gravitinos produced in the same 
process. 

Finally, we shall notice that in models of running mass infiation even the infiaton 
condensate could fragment into non-topological solitons ||8^. In that case Q-balls 
would naturally dominate the Universe. We do not discuss such possibility here. 



6. Summary and Conclusions 

In this paper we have investigated cosmological consequences of models with super- 
heavy gravitinos and/or sfermions. A priori there is no fundamental reason which 
fixes the scale of SUSY breaking. Models with weak-scale SUSY breaking in the 
observable sector have the promise to solve the hierarchy problem. However this 
may not necessarily be the case and the SUSY breaking scale can turn out to be very 
high. Under general circumstances, arbitrarily heavy gravitino mass ^3/2 and/or 
sfermion masses m are quite plausible. Therefore, inspired from the recent models 
of large scale SUSY breaking, it becomes pertinent to re-examine the cosmological 
and phenomenological consequences. 

Gravitino are produced through various processes in the early Universe. Scatter- 
ings of gauge and gaugino quanta in thermal bath, sfermion decays and the infiaton 
decay are the main sources for gravitino production. The main results are presented 
in Eqs. (|2.4|) , (|2.6| ) and ( p.8|) . Sfermion decays usually dominate when the reheating 



temperature Tr > ^ > ^3/2- On the other hand, the contribution from the infiaton 
decay dominates when Tr <^ m. 
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Gravitinos which are heavier than 50 TeV decay before primordial nucleosynthe- 
sis, and hence are not subject to BBN bounds. However, each gravitino produces 
one LSP upon its decay. Hence, in models with conserved i?-parity, the abundance 
of supermassive gravitinos is constrained by the dark matter limit. Indeed, efficient 
annihilation of LSPs produced in gravitino decay sets a lower bound on 1713/2. When 
this lower bound is saturated, gravitino decay can successfully produce non-thermal 
LSP dark matter. 

This is also valid in a gravitino-dominated Universe, which happens when grav- 
itinos are produced very abundantly. However, for 171^/2 > 50 TeV, gravitino dom- 
ination cannot rescue a scenario where thermal LSP abundance at the freeze-out 
exceeds the dark matter bound. For a Wino- or Higgsino-like LSP this is the case 
when the LSP mass > 2 TeV. The reason is that in this case gravitino de- 
cay, while diluting the thermal abundance, leads to non-thermal overproduction of 
LSPs. Therefore, if i?-parity is conserved, gravitinos should never dominate in mod- 
els with such heavy LSPs. The results for gravitino production in conjunction with 
the constraints from the dark matter bound and LSP annihilation are summarized 
in Eqs. (|310|) , (|315|) and {^J^). Figs. (1) and (2) depict the acceptable parts 

of the m — 771^/2 plane for successful scenarios of gravitino non-domination and dom- 
ination respectively. 

We discussed some specific scenarios of baryogenesis in the presence of supermas- 
sive gravitinos. The parameter space for thermal leptogenesis is substantially relaxed 
when 771^/2 > 50 TeV, as a considerably larger reheating temperature Tr and/or right- 
handed (s)neutrino mass Mi will be allowed. This, however, implies that gravitinos 
can also be efficiently produced for a wide range of sfermion masses. Since gravitino 
decay takes place after the completion of leptogenesis, unless they are extremely 
heavy, the generated baryon asymmetry will be diluted in a gravitino-dominated Uni- 
verse. Successful leptogenesis then requires (much) larger right-handed (s)neutrino 
masses than usual. However, it is known that thermal leptogenesis fails for Mi > 10^^ 
GeV since lepton number violating scatterings in this case erase the generated asym- 
metry. This leads to additional constraints on the m — 7713/2 parameter space in 



case of gravitino domination. Our results are summarized in Eqs. ( [4.14| ), ([4.15|) and 
Fig. (4). 

We also considered late time baryogenesis from supersymmetric flat directions 
via the Affleck-Dine mechanism. Thermal effects which can trigger early oscillations 
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of the flat direction condensate, thus suppressing the generate asymmetry, tend to be 
less important for m ^ 1 TeV. A large expectation value for the condensate at the 
onset of its oscillations usually leads to a baryon asymmetry (ne/s) ^ 10~^°, as well 
as a large reheating temperature Tr > m. Gravitinos produced from sfermion decays 
can then dominate the Universe and dilute the initially large asymmetry down to 
acceptable values. The main results in Eqs.( ^.H"D , ( |5.12|) are depicted in Figs (5), (6). 

There is even a closer connection between large baryon asymmetry and efficient 
gravitino production when oscillations of the fiat direction condensate fragment into 
Q-balls (as happen in many cases). Q-balls decay (much) later than the homoge- 
neous condensate, and the larger the baryonic/leptonic charge they carry the longer 
their decay lifetime. Hence the decay of large Q-balls is a natural source for copious 
production of gravitinos which can even dominate the Universe and dilute the large 
baryon asymmetry released by Q-balls. We explicitly demonstrated this for a poten- 
tial with logarithmic corrections, Eqs. ( |5.18|) and ( |5.19|) , but the same conclusions 
hold for other types of fiat potentials. 

To conclude, models with superheavy gravitinos and/or sfermions have very 
interesting cosmological consequences. These models can naturally give rise to a 
large gravitino abundance in the early Universe. This, contrary to models with a 
wealc scale gravitino mass, can turn to a virtue and lead to successful production of 
darlc matter and baryon asymmetry generation. 
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